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Abstract. In this paper, we study cellular automata on Cayley tree of order 2 over the 
(~*) field Z p (the set of prime numbers modulo p). We construct the rule matrix corresponding 

to finite cellular automata on Cayley tree. Further, we analyze the reversibility problem 
of this cellular automata for some given values of a, b, c, d € Z p \ {0} and the levels n of 
Cayley tree. We compute the measure-theoretical entropy of the cellular automata which 



o 



(N 



we define on Cayley tree. 



1. Introduction 



in 

A cellular automaton (plural cellular automata, shortly CA) has been studied and applied 
as a discrete model in many areas of science. Cellular automata (CAs) have very rich 
computational properties and provide different models in computation. CAs were first 
used for modeling various physical and biological processes and especially in computer 
science. Recently, CAs have been widely investigated in many disciplines with different 
\Q ■ purposes such as simulation of natural phenomena, pseudo-random number generation, 

image processing, analysis of universal model of computations, coding theory, cryptography, 
ergodic theory ([D El El HI El E])- 

Most of the studies and applications for CA is extensively done for one-dimensional (1-D) 
CA. "The Game of Life" developed by John H. Conway in the 1960s is an example of a 
two-dimensional (2-D) CA. John von Neumann in the late 40's and early 50's studied CA 
as a self- reproducing simple organisms [7]. 2-D CA with von Neumann neighborhood has 
found many applications and been explored in the literature j8] . Nowadays, 2-D CAs have 
attracted much of the interest. Some basic and precise mathematical models using matrix 
algebra built on field Z2 were reported for characterizing the behavior of two-dimensional 
nearest neighborhood linear CAs with null or periodic boundary conditions [21 El HI El IE] . 
The reversibility problem of some special classes of 1-D CAs reflective and periodic boundary 
conditions has been studied with the help of matrix algebra approach by several researchers 

pun]. 

In Ref. (XT], Fici and Fiorenzi have a first attempt to study topological properties of CA 
on the full tree shift , where Yl* is the free monoid of finite rank | ^ | . In this case the 
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Cayley graph of ^2* is a regular | ^2 |-ary rooted tree. Fici and Fiorenzi [TT] have studied 
cellular automata defined on the full k-axy tree shift (for k > 2). In this paper, we study 
cellular automata on regular Cayley tree of order 2. 

Several notions of entropy of measure-preserving transformation on probability space in 
ergodic theory have been investigated. The notion of entropy, both topological and measure- 
theoretical is one of the fundamental invariants in ergodic theory. In the last years, a lot 
of works have been devoted to this subject [H [12j [131 El- Recall that by the Variational 
Principle the topological entropy is the supremum of the entropies of invariant measures. 
In pQ, the author has shown that the uniform Bernoulli measure is a measure of maximal 
entropy for some 1-D LCAs. Morris and Ward [15] proved that an ergodic additive CA in 
two dimensions has infinite topological entropy (see [16] for details). Recently, Blanchard 
and Tisseur [T7] have introduced the entropy rate of multidimensional CAs and proved 
several results that show that entropy rate of 2-D CA preserve similar properties of the 
entropy of 1-D CA. 

In this short work, firstly we define cellular automata on Cayley tree (or Bethe lattice) of 
order 2. This generalizes the case of one-sided CA (where order of the Cayley tree is one). 
We construct the rule matrix corresponding to finite cellular automata on Cayley tree by 
using matrix algebra built on the field Z p (the set of prime numbers modulo p). Further, 
we discuss the reversibility problem of this cellular automata. Lastly, we study the measure 
theoretical entropy of the CAs on Cayley tree. We show that for CAs on Cayley tree the 
measure entropy with respect to uniform Bernoulli measure is infinity. 

2. Finite CA over Cayley tree 

Let Z p = {0, 1, . . . , p — 1} (p > 2) be the field of the prime numbers modulo p (Z p is 
called a state space). The Cayley tree T k of order k > 1 is an infinite tree, i.e., a graph 
without cycles, from each vertex of which exactly k + 1 edges issue. Let T k = (V, L, i), where 
V is the set of vertices of T k , L is the set of edges of T k and i is the incidence function 
associating each edge lei with its end points x, y G V. A configuration a on V is defined 
as a function x G V — > a(x) G Z p ; in a similar manner one defines configurations a n and 
oj on V n and W n , respectively. The set of all configurations on V (resp. V n , W n ) coincides 
with Q = (resp. tt Vn = Zjf", tt Wn = Zjf"). One can see that tt Vn = x tt Wn . 

Denote by Z^ 2 , i.e., the set of all configurations on V 2 . 

In the sequel we will consider Cayley tree Y 2 = (V, L, i) with the root Xq. If i(£) = {x, y}, 
then x and y are called the nearest neighboring vertices and we write £ =< x,y >. For 
x, y G V, the distance d(x, y) on Cayley tree is defined by the formula 

d(x, y) = min{d\x — x , x\, x 2 , ■ ■ ■ , x^-i, x d = y eV such that the pairs 
< x , X\ >,..., < Xd-i,x d > are neighboring vertices}. 
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W n = {x e V : d(x°, x) = n} 

V n = {x eV : d(x°, x) < n}, 

L n = {i =< x,y >E L : x,y E V n }. 

In this section, we will order the elements of V n in the lexicographical meaning (see [18] ) 
as the Fig. [IJ Given two vertices x u ,x v , the lexicographical order of x u ,x v is defined as 
x u x v if and only if u -< v. 

Let us rewrite the elements of W n in the following order, 



W^ X W n i ■ ■ ■ -> x w n 



One can easily compute equations \ W n \ = 3.2^ _1 ^ and \V n \ = 1 + 3.(2™ — 1). For the sake 
of shortness, throughout the paper we are going to represent vertices x$ , Xy$ , . . . , Xy^™^ 
of W n by means of the coordinate system as follows: 

(1) _ (2) _ (3) _ (4) _ 

X W„ ~ ^11. ..11) x W n ~ X H...12) X w — Xll...21,X w — Xn,,,22, 



(\Wn\-3) _ ((\W n \-2)) _ _([W„|-1)) _ ™(|W„|)) _ 

x Wn — #32.. .211; X Wn — Xja...2\2-,X Wn — X 32 ...221 , X Wn — X 32 ...222- 




x k... ij 



Figure 1. a) Cayley tree of order two with levels 3, b) Elements of the 
nearest neighborhoods surround the center x k ...y, k = 1, 2, 3 and i,j = 1, 2. 

In the Fig. dj we show Cayley tree of order two with levels 3 and the nearest neighborhood 
which comprises three cells which surround the center cell Xk...ij- The state xjjf^y of the cell 
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(i, j) at time (t + 1) is defined by the local rule function / : TJt — > Z p as follows: 



(n T (*+i) _ f(Jt) _(*) _(*) _(*) 

\ L ) x k...ij ~ J \ x k...i> x k...ij' x k...ijl' • h k...ij2J 

= ax kl.ijl + bx klij2 + C 4*L + ^fely ( mod P) > 

where a,b,c,de Z p \ {0}, 4L e W n __ 2 , xg.y e Wn-i and 4% i , 4?..i*2 e A; = 1,2,3 
and i,j = 1,2 (see the Fig. d](b)). 

Specifically, for state 4* in the root vertex we can show 

(2) 4* = f( X 1 ) 4: 4 > X 3 ) = a:C l + ^4 + CX 3^ + f ( m °d J?)- 

Function 

(3) T t : Zf -> Zf 

is called a cellular automaton (CA) generated by the rules (Dp) and (J2]). 

If the boundary cells are connected to 0-state, then CA are called Null Boundary CA, 
i.e., V \ W n = {0} for a fixed n. If the same rule is applied to all of the cells in ever 
evaluation, then those CA are called uniform or regular. 

In Sections [3] and HI we consider linear transformations of finite dimensional vectors 
spaces corresponding to these finite linear cellular automata by imposing the null boundary 
condition, which means that the states of cells outside a given ball around the origin are 
fixed to be zero. 

3. Construction of the rule matrix in the finite case 

In this section, we can characterize finite cellular automata with Null boundary condition 
over Cay ley tree of order two over the field Z p . In order to characterize the corresponding 
rule, first we represent finite Cayley tree n level as a column vector of size (1 + 3(2™ — 1)) x 1. 

Let us denote all configurations of Cayley tree with levels n by Q n . In order to accomplish 
this goal we define the following map 

$ : Q n M (1+3(2 n_ 1))xl (Z p ), 

which takes the t th state X® given by 



O _> yO — (r {t) r {t) r {t) r {t) r {t) r W r W r (t) V 

J£ n t .— IXq , x 1 , . . . , ^21. ..11' x 21...12> • • • ' • L 32...211' x 32...212> x 32...221> x 32...222 ) 



where the superscript T denotes the transpose and M( 1+3 ( 2 ™_i))xi(Zp) is the set of matrices 
with entries Z p . 

The configuration er„ £ is called the configuration matrix (or information matrix) of 
the finite CA on Cayley tree with levels n at time t and is initial information matrix 
of the finite CA. The whole evolution of a particular cellular automata can be comprised 
in its global transition function [5] (see [HI CE21 ED] for the square lattice Z 2 and see [21] for 
the hexagonal lattice). 
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Therefore, one can conclude that $(<7n' ) ) = ^n^ +S (z n -i))xv Using the identification 



due to linearity of the finite CA we can define as follows: 



( M iT ) )(l+3(2"-l))x(l+3(2™-l))^+3 ( 2n_i ))x i 

where n is the number of levels of the Cayley tree. 

Theorem 3.1. Let a, b,c,deZ* p = Z p \{0}, n > 2. Then, the rule matrix {M^ h 
corresponding to the finite cellular automata on Cayley tree of order two with n-level finite 
over NB is given by 

( d P 0i x6 ■■• \ 

Q D 3x3 B 3xe ••• 

C% x3 Dq x q -06x12 ■•• 



X 



(t+1) 



(l+3(2"-l))xl' 



+3(2 n -l)xl+3(2™-l) 
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In the theorem l3.lt we have obtained a general form of the matrix representation (or rule 
matrix) for these linear transformations with respect to a basis given by the lexicographical 
order on the vertices. We do not include a detailed proof of the theorem which gives the 
rule matrix of CA. The proof is obtained by determining the image of the basis elements 
of the space 

z (i+3(2»-i)) under the CA ^ Thege 

images contribute to the columns of the rule 

matrix. 

Let us illustrate this in Examples 13.21 and 13.31 



Example 3.2. If we take the number of level as n = 2, then we get the rule matrix Mr of 
order 10. We consider a configuration of number of levels 2 with null boundary: 




FIGURE 2. A configuration a\ of levels 2 and 3 with null boundary on 
Cayley tree of order two. 
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Let us apply the local rules ([T]) and (J2]) on configuration a 2 in the Fig Then, we get 
a new configuration under this transformation which is 



(t+l) 
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Hence, we obtain the rule matrix M R of order 10 as follows: 
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Example 3.3. Let us consider the configuration with levels 3 given in the Fig. If we 
apply the rules ([I]) and ([2]), then we obtain the following rule matrix: 
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To illustrate the behavior of the finite CA on Cayley tree, we can study the image and 
preimage under finite CA of a configuration by means of relating matrix and its inverse 
matrix (see [TT]). 

4. Reversibility of CA on Cayley tree with Null Boundary 

In this section, we characterize finite cellular automata with NBC determined by nearest 
neighbor rule on Cayley tree. To investigate the reversibility of one dimensional CAs, many 
authors [91 [10] have studied invertibility of the rule matrices corresponding to the CAs. For 
finite CA, in order to obtain the inverse of a finite CA many authors [21 IH El 1201 EI] have 
used the rule matrices. Since we already have found the rule matrix corresponding to 
the the finite CA, by using the matrix in (J4j), we can state the following relation between 
the column vectors and the rule matrix Mr: 

= M^X® (mod p). 
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Table 1 . The reversibility of finite CAs for some given a, b,c,d £ Z* and 
the levels n = 2, 3 of Cayley tree of order two. 
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reversibility of finite LA 
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irreversible 
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irreversible 
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reversible 
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7,11,13,19,23,29 


reversible 



If the rule matrix Mr is non-singular, then we have 

X® — (Mjf ) )~ 1 Jf (t+1) (modp). 

Thus, in this paper one of our main aims is to study whether the rule matrix M# in (j4]) 
is invertible or not. It is well known that the finite CA is reversible if and only if its rule 
matrix is non-singular (see (3j HI [5j [20j [21] for details). If the determinant of a matrix 
is not equal to zero, then it is invertible, so the CA on Cayley tree is reversible, otherwise 
it is irreversible. If the CA is not invertible, then one can study "Garden of Eden" for the 
finite CA (see [III EH]). 

By means of Mathematica, we compute the determinant of the rule matrix for some random 
a,b,c,d £ Z* and the levels n of Cayley tree. 

det{M R 2) ) = -d\c{2{a + b) + c) - d 2 ){-{a + b)c + d 2 ) 2 , 

det{M { R 3) ) = -d 8 {{a + b)c-d 2 f{-2{a + b)c + d 2 ) 2 {{a + b)c 2 {a + b + c)-c{3{a + b) + c)d 2 + d' 1 ). 

We have seen that the CAs are reversible for some given values a, b,c,d £ Z* and n, for 
some values the CAs are irreversible. 

In the Table (TJ we examine under what conditions these linear transformations are in- 
vertible, and check invertibility for a list of parameters using computations by means of 
"Mathematica". For example, if we take &sa = b = c = d = l and n = 3, then we can see 
that the CAs are irreversible for prime numbers p < 47. The reversibility of finite CAs on 
Cayley tree of order two is determined for some given values of a, b,c,d £ Z* and the levels 
n of Cayley tree. One can fully characterize reversibility of finite cellular automata with 
NBC determined by nearest neighbor rule on Cayley tree by computing the determinant of 
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n— 1 

hp(T,a)= lim-^(\/T- 



the matrix in the Eq. (jl]). Also, one can study the reversibility of finite CAs via rank of 
the matrix in the Eq. (jl]) (see [21J). 

5. The measure entropy of the CA on Cayley tree 

In this section we study the measure entropy of cellular automata defined by local rules 
in ([1]) and (j2J) on Cayley tree of order two. In order to state our result, we first recall 
necessary definitions. Let (X, B, /i, T) be a measure-theoretical dynamical system. If a = 
{Ai, . . . , A n } and j3 = {Bi, ... , B m } are two measurable partitions of X, then a V j3 = 
{Ai H Bj : % = 1, . . . , n; j = 1, . . . , m} is the partition of X. Also, T~ 1 a is the partition of 
X and T- l a = {T _: %, . . . , T _1 A n } (see [221 [23] for details). 

Definition 5.1. Let a be a measurable partition of X. The quantity 

H M) = -^2»(A) \og/j(A) 

is called the entropy of the partition a. The logarithm is usually taken to the base 2. Let a 
be a partition with finite entropy, then the quantity 

n— 1 

n— toon 1 v 
i=0 

is called the entropy of a with respect to T. The quantity 

(6) h^{T) = sup{/i M (T, a) : a is a partition with H^(a) < oo} 

a 

is called the measure-theoretical entropy of (X,B, fi,T), the entropy of T (with respect to 
fi). 

Let 7r = {7To, 7Ti, . . . , 7t p -i} be a probability vector. Recall that the Bernoulli measure is 
defined as follows: 

^7r(o[«0, ■ ■ , ik}k) = ^h^i • • - TTifc, 

where q[zq, . . . , is a cylinder set (see [221 123] for details). If we take the Bernoulli measure 
as 

(V ■ i ^ 11 1 1 

AMo «o, • • • , 2fc fcj — — -rir, 

then the measure is called uniform Bernoulli measure, i.e., for all % e Z p , ^(o^]) = ~, then 
H-x is the uniform Bernoulli measure on the space Zj . In this paper, we consider uniform 
Bernoulli measure. 

It is clear that due to (a,p) = l,(b,p) = l,(c,p) = 1 and (d,p) = 1, the rules given in 
the Eqs. ([I]) and (J2J) are bipermutative. The following Theorems have been proved: 

Theorem 5.2. [2U ] Any left-permutative (right-permutative) cellular is surjective (see [25J 
for details). 
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Theorem 5.3. [26J If a cellular automaton is surjective then it preserves a uniform Bernoulli 
measure. 



D'amico et al. jH] have proved that for D- dimensional linear CA with D > 2 the 
topological entropy must be or infinity (see [27]). In the one-dimensional case, the measure 
theoretical entropy of the cellular automata is finite [U [26]. In the following theorem, we 
prove that the linear CA on Cayley tree of order two has infinite entropy. Let us choose 
a, b,c,d G Z* such that the cellular automata Tf defined in the Eq. (J3J) is measure-preserving 
function with respect to (w.r.t.) the uniform Bernoulli measure on the space Z T p . Then we 
have the following theorem. 

Theorem 5.4. Let Tf be cellular automata defined by local rules in (P) and (J2J) on Cayley 
tree of order two over the field Z p . Then the measure theoretical entropy of Tf w.r.t. the 
uniform Bernoulli measure on the space Z^ 2 is infinity. 

Proof. From theorems 15.21 and 15 . 3 [ we note that fj, n is a Tj-invariant measure. Let the zero- 
time partition be given as £(0, 1) = {o[0],o [1], . . . ,o [p — 1]}, we put £(-z,z) = V«=-» cr ~" u £! 
where a is the shift map. Since Tf is permutative, one has 

n-l 

\/Ty fc £(0,l) = £(0,1 + 3(2" - 1)). 

fc=0 

From the definition of measure theoretical entropy w.r.t the measure, we get 

n-l 



V(T /; £(0,1)) = hm-iUVT}- fc £(0,l)) 

k=0 

= -lim- V" fhr(A) log (Mn(A) 

<n — i * 



n— YOOfl 

Aeg(0,l+3(2' l -l)) 



lim I p i+3(2»-i) I log _ J 



n^oon pl+3(2»-l) to p l+3(2"-l) 

lim-(l + 3(2 n - l))logp= oo. 

n— ¥oofl 



Therefore, from the Eq. (Ej), one can conclude that h^iTf) = oo. □ 

Remark 5.5. If we choose the probability vector as % = (1, 0, . . . , 0), then (Tf, £(0, 1)) = 
0. 



6. Conclusions 

In this short paper, firstly we have defined cellular automata on Cayley tree of order 2. 
We have constructed the rule matrix corresponding to finite cellular automata on Cayley 
tree by using matrix algebra built on the field Z p (the set of prime numbers modulo p). 
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Further, we have discussed the reversibility problem of this cellular automata. Lastly, we 
have studied the measure theoretical entropy of the cellular automata on Cayley tree. 

To the best knowledge of the author, it is believed that this is the first instance in the 
literature where such a connection is established. Thus, this connection between cellular 
automata and Cayley tree leads to many questions and applications that waits to be ex- 
plored. Also, investigation of CA on Cayley tree with more higher orders will be studied in 
the future works. 
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